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INSTANTANEOUS ROTATIONAL VELOCITY OF MERCURY 


Han-Shou Liu 
ABSTRACT 


The fluctuation in the angular velocity of the present rotation of 
Mercury is investigated. The instantaneous rotational rate in terms 
of orbital mean motion at different positions along Mercury's orbit is 
given. At aphelion the rotational velocity decreases substantially be- 
cause the solar gravitational torque on the two thermal bulges on Mer- 
cury’s surface tends to retard the rotation of Mercury before aphe- 
lion passage. It is found that the difference between the rotational 
periods derived from the motions at perihelion and aphelion is 4.68 
minutes and that the maximum rate of rotation occurs at f / 4 and 

f = 7-^/4 where f is the true anomaly. 


in 


INSTANTANEOUS ROTATIONAL VELOCITY OF MERCURY 


I. INTRODUCTION 

In the previous paper (Liu, 1969) Liu has suggested that the solar thermal, 
tidal and gravitational effect could account for Mercury's rotation being locked 
into a 3:2 resonance with its orbital revolution. Because of this resonance lock 
the instantaneous rotational period depends strongly on the orientation and posi- 
tion of the planet relative to the Sun. For correct interpretation of the dynamic 
behavior of such a lock-in motion, a precise knowledge about the velocity of 
rotation is of special importance. The present paper is devoted to studying the 
fluctuation of Mercury's rotation in terms of orbital mean motion at different 
positions along its orbit. The calculation is confined to the planar motion. 


II. METHODS OF COMPUTATION 

The rotation of Mercury about its polar axis has the form 

$ = f + 4 > (1) 

where f is the true anomaly and 0, denoting the angle between the axis of the 
two thermal bulges with the radius vector, is the solution of the following differ- 
ential equation (Liu and O'Keefe, 1965) 


d 2 0 

~dF 


2 e sin f 
1 + e cos f 



, 3(B - A) 

C(1 + e cos f) 


cos 0 s in 0 = 0 


( 2 ) 


where (B-A)/C is the present value of the difference in Mercury's equatorial 
moments of inertia and e is the orbital eccentricity. At the present time, the 
influence of the solar thermal and tidal effect on Mercury's rotation is in 
equilibrium. Therefore, the use of equation (2) for orientation determination 
is justified. 


To investigate the fluctuations in the angular velocity of rotation, it is nec- 
essary to transform the variable of function $ from the true anomaly f to the 
mean anomaly M. This can be done from Kepler's equation 


M = n(t - t Q ) = E - e sin E 


(3) 


1 


J3f -ha-" 


where t« = t Q corresponds to the time at perihelidhj n = 2'n/p is the mean motion, 
p - orbital period and E is the eccentric anomaly which is related to 


t an — = 


4 = (i^Y 

2 \1 + cj 


72 f 
t an — 

2 


We c an" Solve equations (2), (3) and (4) on a computer by applying the Ruhge- 
Kutta integration algorithm. The computer program performs the algorithms 
for successively updating expansions. The over-all precision is consistent 
with the 10 -^ truncation level. 


III. RESULTS 

Results of d<t>/dt - fi = (3/2 + S) n, where S is the coefficient of fluctuation, 
are geridr ated for (B-A)/C - 5 x 10 -5 . The intial conditions are chosen such 
that the angle between the axis of the two thermal bulges and the radius vector 
at f = 0 is 0.6° and the averaged value of 0 - 3/2n over the period of revolution 
vanishes. For the orbital eccentricity, the value e = 0.206 is adopted. The 
results of the instantaneous rotational velocity in terms of orbital mean motion 
n at different positions along Mercury's orbit are given in Table 1. It is found 
that 

Q - — n > 0 : 

• " 2 


when I f I < 115° and 


when If I •> 115°. 


n n <0 
2 


The minimum rate of rotation occurs at aphelion because the solar gravita- 
tional torque which exerts a couple on the two thermal bulges on Mercury's sur- 
face tends to retard the rotation of Mercury before aphelion passage. During 
perihelion. passage Mercury's rotation Increases substantially under the effect 
of the Sun's attraction. The maximum rate of rotation occurs at f =45° and 
f = 315°. The dynamic fluctuation in the velocity of rotation is illustrated in 
Figure 1. 


IV. CONCLUSION 

The fluctuation in the rotation of Mercury is calculated accurately to 
iCT 7 of orbital mean motion at the 3 :2 resonance state for which the influence 
of the solar thermal and tidal effect on Mercury's rotation is in equilibrium. 


2 


The difference between the rotational periods derived from the motions at 
perihelion and aphelion is 4.68 minutes. The maximum rate of rotation occurs 
at f = 77/4 and f = 7 tt/ 4. In order to test the theory of rotation for Mercury, 
it seems necessary to improve the accuracy of observations. 
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.'"y, Table 1. Instantaneous Rotational .Velocity at 
r ... . . Different Positions Along Mercury’s Orbit 
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Figure 1-Fluctuation of the Rotational Velocity at Different Positions Along Mercury’s Orbit 
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